Realization of an Excited, Strongly-Correlated Quantum Gas Phase 
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Ultracold atomic physics offers myriad possibilities to study strongly correlated many-body systems in lower 
dimensions. Typically, only ground state phases are accessible. Using a tunable quantum gas of bosonic cesium 
atoms, we realize and control in one dimensional geometry a highly excited quantum phase that is stabilized in 
the presence of attractive interactions by maintaining and strengthening quantum correlations across a confine- 
ment-induced resonance. We diagnose the crossover from repulsive to attractive interactions in terms of the 
stiffness and the energy of the system. Our results open up the experimental study of metastable excited many- 
body phases with strong correlations and their dynamical properties. 



In many-body quantum physics the interplay between 
strong interactions and confinement to a low-dimensional ge- 
ometry amplifies the effects of quantum fluctuations and cor- 
relations. A remarkable example in one dimension is the 
Tonks-Girardeau (TG) gas, where bosons with strong repul- 
sive interactions minimize their interaction energy by avoid- 
ing spatial overlap and acquire fermionic properties JT] [2j. 
Evidence for this ground state phase was found using Bose- 
Einstein condensates (BEC) loaded into optical lattices 13] |4). 
While many-body quantum systems are usually found in their 
ground state phases, long-lived excited state phases are re- 
sponsible for some of the most striking physical effects, ex- 
amples ranging from vortex lattices in superfluids to subtle 
topological excitations in spin liquids 0. However, the ex- 
perimental realization of excited phases is difficult, as these 
usually quickly decay by intrinsic effects or by coupling to 
the environment. In this context, cold atoms [3, 4.I6 HT21 may 
provide unique opportunities for the realization of long-lived, 
strongly interacting, excited many-body phases due to the ex- 
cellent decoupling from the environment and the tunability of 
interactions via, for example, Feshbach resonances. 

For an ultracold one-dimensional (ID) system of bosons, 
we prepare a highly-excited many-body phase known as 
the super- Tonks-Girardeau (sTG) gas ff3l . In this highly- 
correlated quantum phase, interactions are attractive, and 
rapid decay into a cluster-type ground state is in principle pos- 
sible. However, a surprising property of this many-body phase 
is its metastability. Attractive interactions strengthen correla- 
tions between particle positions and ensure, similar to an ef- 
fective long-range repulsive interaction, that particles rarely 
come together. To realize this exotic phase, we observe and 
exploit a ID confinement-induced resonance (CIR) lfl4l [131 . 
This resonance allows us to first enter deeply into the repulsive 
TG regime to establish strong particle correlations and then to 
switch interactions from strongly repulsive to strongly attrac- 
tive. The frequency ratio of the two lowest-energy collective 
modes [16| provides accurate diagnostics for the crossover 
from the TG to the sTG regime. In particle loss and expan- 
sion measurements we study the time evolution of the system 



through the crossover. 

We tune the strength of the interaction as characterized by 
the three-dimensional (3D) scattering length a 3 D by means of 
a magnetically-induced Feshbach resonance [17]. For a ID 
system, a CIR arises and strongly modifies the ID scatter- 
ing properties when a 3 rj approaches the harmonic oscillator 
length aj_ = dWj (moj±) of the transversal confinement with 
trap frequency oj± lfl4l [T5l . Here, m is the mass of the parti- 
cles and h is Planck's constant divided by 2ir. More precisely, 
the coupling constant giD of the ID <5-function contact poten- 
tial Uid{z) = g\Y)8{z) behaves as lfl4ll 

2h 2 2h 2 a 3D 1 

9id — 2~ i n / — ' y l ) 

maiu ma± 1 — C a 3 D/a_i_ 

where aiD is the ID scattering length defined by this equation 
and C = 1.0326 is a constant. Thus, the CIR allows tuning of 
<7iX). For values of a 3 D less but close to a±/C (a 3 D %ax/C) 
the coupling parameter guy is large and positive, and for 
a 3 D > aj_ jC it is large and negative, leading to an effectively 
attractive interaction. For homogenous systems with gio > 0, 
it is customary to characterize the strength of interactions by 
the Lieb-Liniger parameter 7 = giT>m/(H 2 rirD), where niD 
is the linear ID density of the system ||2l|6). The TG gas cor- 
responds to the limit 7 > 1 01 — > 00. As interactions are 
increased, the system becomes strongly correlated and is fully 
dominated by its kinetic energy. In previous experiments, 
without the capability to tune a 3 D, a maximum of 7 rj 5.5 
was achieved [4|, while an effective strength 7 ff s» 200 was 
reached with an additional shallow lattice potential along the 
longitudinal direction [3 |. In the former experiment, a satura- 
tion for the size and energy of the ID system was observed, 
and in the latter experiment the momentum distribution was 
studied. 

But what happens in the case of strong attractive interac- 
tions <7id — > —00, i.e. aiD > 0? The ground state for a 
system of N attractively interacting bosons in ID is a cluster 
state lfl8l [T9l . which one would expect, in a cold atom sys- 
tem, to decay quickly via molecular channels. However, by 
crossing the CIR from the TG side, i.e. switching interac- 
tions from gin — +00 to gm — —00, an excited gas-like 
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FIG. 1: A, Experimental setup. The lattice potential is created by 
two retro-reflected laser beams confining the atoms to an array of 
one-dimensional tubes with equipotential surfaces shown in red. B, 
Along each tube (left) we excite the lowest compressional mode 
(center) and compare its frequency to the dipole mode (right). C, 
The strength of the interatomic interaction is adjusted by tuning the 
s-wave scattering length a^. The background scattering length rises 
gently from to 1240 ao when the magnetic field B is tuned from 
17 to 76 G. Further tuning is possible near a Feshbach resonance at 
47.78(1) G to absolute values beyond 4000 ao. The dashed line indi- 
cates a±/C for a transversal trap frequency of uj± — 2n x 13. 1 kHz. 
D and E present typical data sets for the compressional mode in the 
TG and sTG regime at a 3D = 875(1) a and a 3D = 2300(200) a , 
respectively. The upper panels show the atom number, the lower 
panels show the 1/e-cloud- width after time-of-flight. The solid 
lines in the lower panels are sinusoidal fits (see online material), 
yielding the oscillation frequencies u)c = 2n x 30.6(3) Hz and 
loc = 27T x 241(1) Hz, respectively. 



phase, the sTG gas, should be accessible lfT3ll . Is this ex- 
cited phase stable, i.e. does it exist at all? The expectation 
is that the large kinetic energy inherited from the TG gas, in a 
Fermi-pressure like manner, prevents the gas from collapsing 
||20l . This stability can most simply be inferred from a Bethe- 
ansatz solution to the Lieb-Liniger model with attractive inter- 
actions |20||2TI . This ansatz yields for the energy per particle 
E/N « h 2 TT 2 n 2 B /[6m(l — niD^m) 2 ], corresponding to the 
energy of a gas of hard rods (T), for which any represents the 
excluded volume. This results in a positive inverse compress- 
ibility and also in an increased stiffness of the systems as long 
as niDaiD is sufficiently small. Interestingly, in this phase the 
density correlations are even stronger than in the TG gas, as 
they show a power-law decay that is slower than for a TG gas 
ifTSl . indicating an effective long-range interaction. 

We realize the crossover all the way from a non-interacting 
gas via the ID mean-field Thomas-Fermi (TF) regime to a 
TG gas and then to a sTG gas. We exploit the fact that our 
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FIG. 2: Transition from the non-interacting regime via the mean- 
field TF regime into the TG regime. The squared frequency ratio 
R = icq/iJjj of the lowest compressional mode with frequency ujc 
and the dipole mode with frequency ujd serves as an indicator for 
the different regimes of interaction. For increasing interactions from 
7 = to 7 w 500 the system passes from the ideal gas regime (R = 
4) to the ID TF regime (R « 3) and then deeply into the TG regime 
(R = 4). The inset shows the transition from the non-interacting 
regime to the mean-field regime in more detail. The vertical error 
bars refer to standard error and the horizontal error bars reflect the 
uncertainty in determining am and nm (see online material). The 
horizontal error bar on the data point at 7 = (not shown in the 
inset) is ±0.03 ao. 



ID systems possess weak harmonic confinement along the 
axial direction characterized by the confinement length an. 
Whereas the frequency uijof the lowest dipole mode depends 
only on the confinement, the frequency ljq of the lowest ax- 
ial compressional mode is sensitive to the various regimes of 
interaction [16|. For the non-interacting system one expects 
R = ujq/uj 2 -) = 4. This value then changes to R — 3 for 
weakly repulsive interactions in a ID TF regime [7|. For in- 
creasing positive interaction strength, R is expected to change 
smoothly to 4 when entering the TG regime as the system be- 
comes fermionized and hence effectively non-interacting. A 
rise beyond the value of 4, after crossing the CIR, would then 
constitute clear evidence for the sTG regime [ 13 1. As ayo is 
further increased, the system will finally become unstable and 
R is expected to turn over and drop towards zero. For a har- 
monically confined system, the point of instability is reached 
when the overall length of the system of hard rods, Nam, be- 
comes of the order of the size \/N ay for the wave function of 
N non-interacting fermions, i.e. A = N a\j) / (\N~ a\\) ~ 1. 
We use A 2 as an alternative parameter to 7 to characterize 
the strength of the interaction as it accounts for the harmonic 
confinement. 



We start from a 3D Bose-Einstein condensate (BEC) with 
up to 2 x 10 5 Cs atoms with no detectable thermal fraction 
in a crossed-beam dipole trap with magnetic levitation 11221 . 
Depending on the interaction regime to be studied, we then 
set the number of atoms in the BEC to values in the range of 
(1 — 4) x 10 4 by means of forced radio-frequency evapora- 
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FIG. 3: A, The ratio R — loq/uj 2 ^ is plotted as a function of the 
interaction parameter A 2 — Na 2 D / a 2 . The squares show the mea- 
surements in the attractive regime (giD < 0), providing evidence 
for the super-Tonks-Girardeau gas. The circles show the transition 
from the TF to the TG regime (<?id > 0, same data as in Fig. 2 
for 7 > 1). The solid (dashed) line presents the theoretical data for 
<7id > (<?id < 0) by Astrakharchik et al.|13|. The dotted line 
corresponds to the model of hard rods. For reference, the measure- 
ments for (/id < are numbered. Data points lc to 6 are taken at 
ojd = 2tv x 115.6(3) Hz. For data points la and lb the trap fre- 
quency is ujd = 27r x 22.4(1) Hz and ujd = 2n x 52.3(1) Hz, re- 
spectively. For all measurements in the sTG regime a± = 1346(5) 
ao. B, The parameters azu (dashed-dotted), aro (solid), and gw 
(dashed) are plotted in the vicinity of the Feshbach resonance (FR) at 
47.78(1) G. The horizontal dotted line indicates the value of a±/C. 
The pole of the CIR is at 47.36(2) G. a 1D has a pole (P) at 47.96(2) 
G. The bell-shaped curve at the bottom left indicates the atomic dis- 
tribution as a function of the magnetic field determined from high- 
resolution microwave spectroscopy. 



tion. To confine the atoms in ID, i.e. to freeze out transversal 
motion, we use a two-dimensional optical lattice 11211 . which 
forms an array of vertically oriented elongated tubes with an 
aspect ratio that we set to values between 100 and 1000 (Fig. 
1 A). We occupy between 3000 — 6000 independent tubes with 
8-25 atoms in the center tube. The interaction strength <?id is 
controlled by magnetic tuning of a^n by means of a combi- 
nation of a broad and a narrow Feshbach resonance (Fig. 1C) 
with poles at B = -11.1(6) G and B — 47.78(1) G and 



widths of about 29.2 G and 164 mG, respectively ll23l . The 
broad resonance provides a slow variation of a 3 D, allowing us 
to gently tune 030 from ao near 17.119 G to about 1240 ao 
near 76 G, while the narrow resonance allows us to tune a 3 o 
to absolute values beyond 4000 ao given our magnetic field 
control. We convert the applied magnetic field B into us- 
ing the fit formula of Ref. 23. A magnetic field gradient, used 
to levitate the atomic sample [24|, introduces a small spread 
in the value of across the sample. 

To determine the oscillation frequencies u>c and ujd of the 
fundamental modes (Fig. 1 B), we excite each mode sepa- 
rately at a given value of the magnetic field B [24] and let 
the atoms evolve for a varying amount of hold time. The dis- 
tribution is then imaged in momentum space by taking an ab- 
sorption picture after release and expansion. To avoid possible 
broadening effects due to interaction during the initial expan- 
sion, CI3D is set to zero near B = 17.119 G at the moment 
of release. To extract the frequency, we determine for each 
hold time the axial 1 /e-width of the distribution and then fit a 
damped sinusoid with linear offset to this data. Typical mea- 
surements of u>c are shown in Fig. 1 D and E. Whereas the 
atom number remains constant for <?id > 0, we observe some 
atom loss and a slight broadening of the distribution for attrac- 
tive ID interactions. In all parameter regimes, the ID system 
is sufficiently stable to allow a reliable measurement of u>c- 

First, we show that we can tune the system from the non- 
interacting regime deeply into the repulsive TG regime (Fig. 
2). In agreement with expectations, the value for R = ojfj/ujfj 
first drops from 4 to 3 and then increases back to 4 as 7 
is tuned by means of the gently-varying background scatter- 
ing length. We find that the TG regime is fully reached for 
7 > 50. A further increase to values up to 7 m 500 does 
not lead to changes for R. Note that, as (I3D approaches a±, 
the divergence of gm according to Eq. 1 has to be taken into 
account when determining 7 11241 . Heating of the system can 
be excluded as we can return to a 3D BEC without significant 
thermal background when ramping down the lattice potential. 

The attractive regime is entered by crossing the CIR on the 
low-field wing of the 47.78 G Feshbach resonance, oid is now 
small and positive. The central results of this work are sum- 
marized in Fig. 3A and compared to the theoretical work of 
Ref. 13. We plot R = uj"q/u: 2 d as a function of the interaction 
parameter A 2 . For reference, Fig. 3B plots a 3 o, aiD, and g^ 
in the vicinity of the Feshbach resonance as a function of the 
magnetic field B. As the CIR is crossed and A 2 is increased, 
R rises beyond the value of 4. This provides clear evidence for 
the sTG regime as R = 4 is the maximal value for bosons with 
repulsive contact interaction. This increase is expected from 
the model of a gas of hard rods, and our data initially follows 
the prediction from this model. However, as A 2 is increased, 
R reaches a maximum and then starts to drop. The maximum 
of about 4.5 is reached for A 2 s» 3 x 10 -2 . The existence 
of the maximum is in qualitative agreement with the results 
obtained from Monte-Carlo simulations [13]. The theoretical 
prediction, however, underestimates the measured R. This is 
probably due to the local density approximation, which may 
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FIG. 4: Stability and kinetic energy in the TG and sTG regimes. A, 
relative number of atoms remaining and B, relative 1/e- width along 
the axial direction after 10 ms expansion, after a hold time r = 10, 
50, 100, and 200 ms (circles, triangles, squares, and diamonds, re- 
spectively) at a given magnetic field B. The position of the CIR, 
the pole of the Feshbach resonance (FR), and the pole for am (P) 
are as indicated. For these measurements a± = 1523(6) ao and 
ujd = 2-7T x 115.6(3) Hz. The atom number is normalized to the 
initial value of 1.7(1) x 10 4 and the width is normalized to the initial 
value in the TG regime. 



not be applicable to our system with low particle numbers. For 
comparison, the results from Fig. 2 for 7 > 1 are shown. Note 
that 7 w 500 corresponds to small values of A 2 « 10 -4 . For 
this data, at higher particle numbers, there is excellent agree- 
ment with the theoretical prediction (solid line) in the entire 
crossover from the mean-field regime to the TG regime] 16|. 

We study the stability of the system in the crossover from 
the TG to the sTG regime and find further evidence for the 
existence of the CIR by recording particle loss and measur- 
ing the axial width of the atomic cloud after release from the 
tubes. The axial width is a measure for the kinetic energy of 
the system as interactions are instantly switched off upon re- 
lease. Similar conditions are used as for the measurements on 
the sTG regime presented in Fig. 3. The TG regime is entered 
adiabatically to avoid the excitation of collective modes. The 
system is prepared at a^u — 887(1) ao at a magnetic field of 
B = 42.77(2) G with about 11 atoms in the central tube. The 
magnetic field is then ramped to a specific value within 0.2 
ms and the sample is held at this value for a variable hold time 
r from 10 to 200 ms. a± is set to 1523(6) ao. The results 



(Fig. 4) for different hold times r in the tubes show that, for 
t = 10 ms, corresponding to the timescale of the measure- 
ments in the sTG regime shown in Fig. 3, the transition from 
the TG to the sTG regime appears very smooth. There is es- 
sentially no particle loss when the system is deep in the TG 
regime and close to the CIR. The loss gradually increases in 
the attractive regime as one moves to larger values of B and 
towards the pole for am. Correspondingly, the width of the 
sample exhibits a smooth behavior across the CIR, showing a 
slight increase for larger B. This behavior is consistent with 
the expectation of an increased energy in the sTG regime lfl3l . 

For longer hold times, the data for the atom number and the 
sample width develop distinct features at the calculated posi- 
tion of the CIR. Evidently, the system is in a transient state. 
For r = 50 ms, the number of remaining atoms shows a dip 
that correlates with a peak in the kinetic energy of the sample. 
Both features become more prominent and asymmetric for 
longer hold times (t = 100 and 200 ms). Note that, in com- 
parison, no pronounced effects are visible at the pole of the 
Feshbach resonance for a^u ■ Our results must be connected 
to the fact that the energy spectrum of the system changes dra- 
matically across the CIR, from the TG to the sTG regime fl9l . 
The system acquires a deeply lying ground state together with 
a family of lower lying many-body excited states, potentially 
opening up new decay channels. Also, the CIR strongly mod- 
ifies the two-body scattering problem, making formation of 
confinement-induced molecules in trans vers ally excited trap 
states [14] possible. 

The non-trivial time evolution observed in this system 
raises intriguing questions on possible coupling and decay 
mechanisms for strongly interacting excited many-body sys- 
tems, in particular in the context of integrability of ID systems 
||25ll . Our results offer an example of the counter-intuitive ef- 
fects that occur in many-body systems, and open up the possi- 
bility to study the dynamical properties of strongly-correlated 
systems with effective long-range interactions ll26l l27ll un- 
der conditions where all parameters are tunable and, in fact, 
can be changed dynamically. Similar to magnetic Feshbach 
resonances in atomic scattering, we expect the confinement- 
induced resonance demonstrated here to serve as a general 
tool to tailor interactions in ID and possibly also in 2D sys- 
tems 11281 . allowing for the further investigation of strongly 
correlated phases in the context of cold atomic gases. 
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Materials and methods 



center of the array. 



Lattice loading. We produce a BEC of Cs atoms in the 
lowest hyperfine sublevel with hyperfine quantum numbers 
F = 3 and m F = 3 in a crossed beam dipole trap with trap 
frequencies oo x , VtZ = 2n x (15,20, 13) Hz, where z denotes 
the vertical direction. The BEC is adiabatically transferred 
from the dipole trap to the array of tubes by exponentially 
ramping up the power in the lattice laser beams with waists 
<~ 350 /im within 500 ms. The repulsive interaction causes 
the atoms to move radially outwards during the initial phase of 
the lattice loading in response to the strong local compression. 
We use this effect to vary the total number of tubes loaded and 
hence the atom number per tube by setting a 3 D for the load- 
ing process to values between 40 ao and 350 ao- For the data 
set in the repulsive regime (Fig.3A, circles), we exponentially 
ramp down the crossed beam dipole trap during the loading 
process and reach longitudinal and transversal trap frequen- 
cies of uj d = 2-7T x 15.4(1) Hz and loj_ = 2ir x 13.1(1) kHz 
with a transversal confinement length a± = 1440(6) a . Here, 
depending on the regime of interaction to be studied, the num- 
ber of atoms in the central tube is set to values between 8 and 
25. For the data set in the sTG regime (Fig. 3 A, squares) we in- 
crease cod to 2n x 115.6(3) Hz to reduce the vertical extent of 
the sample and hence the variation of the magnetic field across 
the atom cloud, see below. For this, we keep the depth of the 
crossed beam dipole trap constant during the loading process 
and then ramp up the power in one of the beams within 100 
ms. In this regime we choose uj± — 2ir x 15.0(1) kHz, cor- 
responding to a± = 1346(5) ao- The number of atoms in the 
central tube is set to values between 8 and 1 1 . 

Array of ID tubes. The atom number per tube becomes 
fixed once tunneling is suppressed during the loading process 
and can be determined by integrating a Thomas-Fermi profile 
along the tubes (1). The number of atoms in tube is 
given by 



N, 



N , 



.dlat 



.d\ a t 



1 3/2 



N ,o 



2lT R x Ry 



where AT to t is the total atom number, No,o is the occupation 
of the central tube, d lat = A/2 is the lattice spacing at a wave- 
length A = 1064.5 nm, and R x ^ y are the Thomas-Fermi radii 
in the horizontal directions. To calculate the effective atom 
number per tube N, we average over the tubes, weighting each 
tube by its atom number. This procedure accounts for the fact 
that we measure an averaged frequency u>c, as loq is expected 
to slightly vary from tube to tube. The result for lo c should 
be dominated by the more heavily occupied tubes close to the 



Magnetic levitation. To hold the (F = 3, mp = 3) atoms 
in the vertically oriented tubes, magnetic levitation by means 
of a magnetic field gradient of 3 1 . 1 G/cm is applied. The gra- 
dient introduces a small field spread over the atomic sample. 
This sets our precision to tune the interaction strength. For 
the measurements in the sTG regime the distribution has a full 
width at half maximum (FWHM) of 30 mG. We measure the 
atom distribution in the magnetic field by driving a magnetic- 
field-dependent microwave transition. A typical distribution 
is shown in the bottom left corner of Fig.3B. 

Excitation of collective modes. We use two different 
methods to excite the lowest compressional mode. For a mea- 
surement in the mean-field regime, we use a rapid change of 
the interaction strength to excite the oscillation. For this, we 
ramp the scattering length adiabatically in 100 ms to a value 
that is approximately 50 a from the desired final value and 
then perform the last part of the ramp non-adiabatically. For a 
measurement in the sTG regime, we use an analogous method. 
We simply ramp sufficiently quickly, within about 5 ms, all the 
way from the mean-field across the TG into the sTG regime. 
For the TG regime, we chose to excite the mode by compress- 
ing the cloud adiabatically with an additional dipole trap laser 
beam, starting the motion by rapidly ramping down the power 
of this beam. In all cases, we adjust the ramp speeds so that the 
measured oscillation amplitude is within 10-20% of the initial 
cloud size. To excite the dipole mode at frequency cod, we 
adiabatically lower the levitating magnetic field gradient and 
hence displace the cloud along the vertical direction. Quickly 
readjusting the gradient back to full levitation leads to excita- 
tion of the dipole oscillation. 

Determination of 7. We make a conservative estimate to 
determine the Lieb-Liniger interaction parameter 7 



7 = 
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To take into account that the atom number varies according to 
Nij, we first calculate for every tube separately. We cal- 
culate the center density for each tube both in the mean-field 
and in the TG regime and use the larger value to determine 
7i j. We then take 7 as the weighted average over 7^. The 
error in determining 7 largely comes from the determination 
of din, reflecting the magnetic field distribution across the 
sample. 
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